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o 

, Diabatic description of rotational bands provides a clear-cut picture for understanding 

the back-bending phenomena, where the internal structure of the yrast band changes dramat- 
ically as a function of angular momentum. A microscopic framework to obtain the diabatic 
bands within the mean-field approximation is presented by making use of the selfconsistent 
collective coordinate method. Applying the framework, both the ground state rotational 
bands and the Stockholm bands are studied systematically for the rare-earth deformed nu- 
clei. An overall agreement has been achieved between the calculated and observed rotational 
spectra. It is also shown that the inclusion of the double-stretched quadrupole-pairing in- 
^ . teraction is crucial to obtain an overall agreement for the even-odd mass differences and the 

j^-^ ' rotational spectra simultaneously. 

o 

^ . §1. Introduction 

o 

' Back-bending of the yrast rotational bands is one ofthe most striking phenomena 

in the spectroscopic studies of rapidly rotating nuclei. HP 'H-* The first back-bending, 
which has been observed systematically in the rotational bands of the rare-earth 
nuclei, has been understood as a band-crossing between the ground state rotational 
^ . band {g-hand) and the lowest two quasineutron excited band (s-band). A simple 

\ approach to describe the band-crossing is the cranked mean-field approximation, 

where the concept of independent particle motion in the rotating frame is fully 
employed. As long as the conventional (adiabatic) cranking model is used, however, 
. , the two bands mix at the same rotational frequency and, in the crossing region, loose 

' their identities as individual rotational bands. It should be noted that the difficulty 

lies in the fact that the angular momenta of two bands are considerably different in 
the vicinity of the crossing frequency where the mixing takes place, esDeciaUy in the 
case of sharp back-b endings, and such a mixing is largely unphysical.U'*'Q)'B) 

A key to solve this problem is to construct diabatic rotational hands, where the 
internal structure of the band does not change abruptly. %%%% Once reliable 
diabatic bands are obtained it is rather straightforward to mix them if the number 
of independent bands are few as in the case of the first back-bending. Note, however, 
that it is highly non-trivial how to construct reliable diabatic bands in the mean-field 
approximation, because it is based on the variational principle and the mixing at the 
same rotational frequency is inevitable for states with the same quantum numbers (in 
the intrinsic frame). On the other hand, in the mean-field approximation, the effects 
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of rotational motion on the internal structure of the g-hand can be nicely taken into 
account as selfconsistent changes of the deformation and the pairing gap parameters. 
Furthermore, rotation alignment effects of the quasiparticle angular momenta are 
described in a simple and clear way. Therefore, it is desired to develop a method to 
describe the rotational band diabatically within the mean-field approximation. 

In this paper, we present a powerful method to obtain reliable diabatic rotational 
bands by making use of the selfconsistent collective coordinate (SCC) method. 
The method is applied to the g- and s-bands and the results for nuclei in the rare- 
earth region are compared systematically with experimental data. In order to repro- 
duce the rotational spectra, the choice of residiial interaction is essential. We use the 
pairing-plus-quadrupole force type interaction. HIP However, it has been well-known 
that the moment of inertia is generally underestimated by about 20—30% if only 
the monopole-pairing interaction is included. Ii3 Therefore, we exploit the monopole 
and quadrupole type interaction in the pairing channel, and investigate the best 
form of the quadrupole-pairing part. This is done in §2. After fixing the suitable 
residual interaction, we present in §3 a formulation to describe the diabtic rotational 
bands and results of its application to nuclei in the rare-earth region. In practical 
applications it often happens that a complete set of the diabatic quasiparticle basis is 
necessary; for example, in order to go beyond the mean-field approximation. For this 
purpose, we present in §4 a practical method to construct the diabatic quasiparticle 
basis satisfying the orthonormality condition. Concluding remarks are given in §5. 

§2. Quadrupole-pairing interaction suitable for deformed nuclei 

In this section we try to fix the form of residual interactions, which is suitable 
to describe the properties of deformed rotating nuclei. It might be desirable to use 
effective interactions like Skyrme-type interactions, B but that is out of scope of the 
present investigation. We assume the separable-type schematic interactions instead, 
and try to fix their forms and strengths by a global fit of the basic properties; the 
even-odd mass difference and the moment of inertia. 

2.1. Residual interactions 

The residual interaction we use in the present work is of the following form: 

V = -GoPjo^oo -G2J2 pIkP2K -\Y1 ^2kQ\kQ2K. (2-1) 

K K 

where the first and the second terms are the monopole- and quadrupole-pairing in- 
teractions, while the third term is the quadrupole particle-hole type interaction. The 
pairing interactions are set up for neutrons and protons separately (the T = 1 and 
= ±1 pairing) as usual, although it is not stated explicitly, and only the isoscalar 
part is considered for the quadrupole interaction. The quadrupole-pairing interac- 
tion is included for pujmose of better description of moment of inertia: It has been 
known for many yearsUB) that the cranking moments of inertia evaluated taking into 
account only the monopole-pairing interaction underestimate the experimental ones 
systematically in the rare-earth region, as long as the monopole-pairing strength 
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is fixed to reproduces the even-odd mass differences. It sliould be mentioned that 
the treatment of residual interactions in the pairing and the particle-hole channels 
are different: In the pairing channel the mean- field (pairing gap) is determined by 
the interaction selfconsistently, while that in the particle-hole channel (spatial de- 
formation) is obtained by the Nilsson-Strutinsky methodEJ^'HM) and the interaction 
in this channel only describes the dynamical effects, i.e. the fluctuations around the 
equilibrium mean-field. 

The basic quantity for deformed nuclei is the equilibrium deformation. For the 
present investigation, where the properties of deformed rotational nuclei are system- 
atically studied, the Nilsson-Strutinsky method is most suitable to determine the 
equilibrium deformations, because there is no adjustable parameters. As empha- 
sized by Kishimoto and Sakamoto, HIP the particle-hole typenimdmpole interaction 
for deformed nuclei should be of the double-stretched form: 113' 'EJ^'tEP 

Q2K = E ^^K{ij) c\c,, q2K{ij) = {i\{r^Y2K)"\j), (2-2) 

where c| is the nucleon creation operator in the Nilsson state (O)" means that 
the Cartesian coordinate in the operator O should be replaced such as ^ X I, = 
{LOk/uJo)xk {k = x,y,z), where lOx, oJy and ujz are frequencies of the_anisotropic 
oscillator potential and related to the deformation parameter (e2, 7);ll^'"^ here 
= 41.0/^1/3 MeV {A is the mass number). Then the selfcon- 
sistent condition gives, at the equilibrium shape, a vanishing mean value for the 
double-stretched quadrupole operator, {Q2k) = 0, and thus the meaning of resid- 
ual interaction is apparent for the double-stretched interaction. Moreover, the force 
strengths are determined at the same time to be the so-called selfconsistent value, 

by which the /?- and 7-vibrational excitations are correctly described. Strictly speak- 
ing, the vanishing mean value of Q2K holds only for the harmonic oscillator model. 
It is, however, easily confirmed that the mean value vanishes in a good approxima- 
tion in the case of Nilsson potential. In fact the calculated ratio of mean values of 
the double-stretched and non-stretched quadrupole operator is typically within few 
percent, if the deformation parameter determined by the Strutinsky procedure is 
used. 

Pairing correlations are important for the nuclear structure problem as well. The 
operators entering in the pairing type residual interactions are of the form 

^00 = E 44' 4x = E P2Kiij) 44 (2-4) 

j>0 ij>0 

where j denotes the time-reversal conjugate of the Nilsson state j. In contrast to 
the residual interactions in the particle-hole channel, there is no such selfconsitency 
condition known in the pairing channel. Therefore, we use the Hartree-Bogoliubov 
(HB) procedure (exchange terms are neglected) for the pairing interactions, only 
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the monopole part of which leads to the ordinary BCS treatment. Note that the 
generalized Bogoliubov transformation is necessary in order to treat the quadrupole- 
pairing interaction, since the pairing potential becomes state-dependent and contains 
non-diagonal elements: 

Aij = Aqo 5ij + ^ A2KP2K{ij), (2-5) 

K 

where Aqq = Gq {Pqo) and A2K = G2 {P2k), the expectation values being taken with 
respect to the resultant HB state. 

For the application of these residual interaction we are mainly concerned with 
deformed nuclei in the rare-earth region, where the neutron and proton numbers 
are considerably different. In such a case, the "iso-str etching" of multipole oper- 
ators, £21)' Qr {2Nr/Afl^Qr for T = v^TT [N-r denotes the neutron or pro- 
ton number and A the mass number), is necessary in accordance with the differ- 
ence of the oscillator frequencies, uoq = {2Nt- / AY/^lvq, or of the oscillator length, 
(6g)^ = {2Nr / A)~^/^h'Q. We employ this modification for the quadrupole interaction 
in the particle-hole channel. 

2.2. Treatment of pairing interactions 

As for the quadrupole-pairius part, there are at least three variants that have 
been used in the literature. cHP'E3)'EJ)'E3)'E3^'E3) Namely, they are non-stretched, 
single-stretched and double-stretched quadrupole-pairing interactions, where the pair- 
ing form factor in the operator in Eq. ( |2-4| ) is defined as, 

V2K{ij) = ii\r^Y2K\j), (i\{r^Y2K)'\3), (i|(r'l2/<)"|j), (2-6) 

respectively. The single-stretching of operators is analogously performed by the 
replacement, Xk — > x'^ = \/uJkh^ {k = x,y,z). Note that there are matrix ele- 
ments between the Nilsson states with zlA'^osc = ±2 in Eq. ( |2-6| ). We have neglected 
them in the generalized Bogoliubov transformation in accordance with the treat- 
ment of the Nilsson potential, which is arranged to have vanishing matrix elements 
of Z\7Vosc = ±20 

Consistently to the Nilsson-Strutinsky method, we use the smoothed pairing gap 
method in which the monopole-pairing force strength is determined for a given 
set of single-particle energies by 



J- = i log [a/ A + ^J{A/AY + l^ , (2-7) 

where g-p is the Strutinsky smoothed single-particle level density at the Fermi surface, 
A is the cutoff energy of pairing model space, for which we use A = 1.2^a;o, and A 
is the smoothed pairing gap. We introduce a parameter d (MeV) to control the 
strength of the monopole-pairing force by 



*■* The hexadecapole deformation leads extra Z\A'^osc = ±2, ±4 coupling terms, but they are 
neglected in our calculations. 
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through Eq. (2-7), where the same smoothed pairing gap is used for both neutrons 
and protons, for simpUcity. As for the quadrupole-pairing force strength, we take 
the following form, 

G2 = Go with Ro = \.2A^I^ fm. (2-9) 

Thus, we have two parameters d (in MeV) and 52 for the residual interactions in the 
pairing channel. 

It is worthwhile mentioning that Eq. (|2-7D gives the form. 



G5-^(^) , with c=^,7^1og(2^/^), (2.10) 

for the semiclassical treatment of the isotropic harmonic oscillator model, 
where ~ (3A^)^/^/(/ia'o)) and it is a good approximation for the Nilsson po- 
tential. The quantity log depends very slowly on the mass number and can 
be replaced by a representative value for a restricted region of mass table. Taking 
A = 170, one obtains c « 23, which gives the monopole-pairing force strength often 
used for nuclei in the rare-earth region. 

2.3. Determination of parameters d and g2 

Now let us determine the form of the quadrupole-pairing interaction. Namely, we 
would like to answer the question of which form factor in Eq. ( |2-6| ) is best, and of what 
are the values of the parameters, d and §2, introduced in the previous subsection. 
For this purpose, we adopt the following criteria; the moments of inertia J7cLi{f the 
Harris formulae^' and the even-odd mass differences (the third order formulaEs*) for 
even-even nuclei should be simultaneously reproduced as good as possible. Since the 
neutron contribution is more important for the moment of inertia, we have used the 
even-odd mass difference for neutrons, E^''^\ Then it turns out that the proton 
even-odd mass difference is also reasonably well reproduced as long as the same 
smoothed pairing gap is used for neutrons and protons. Thus, the two parameters 
d (MeV) and 52 are searched so as to minimize the root-mean-square deviations of 
these quantities divided by their average values. 



1 -'Vdata , , , , 1/2 / r 1 "dut: 



(exp) 



(2-11) 



for X = Jq and e\j'°^ . Nuclei used in the search are chosen from even-even rare-earth 
nuclei in Table |, thus A^data = 83 (58) for x = Ei""'"^ (Jo)- 

Table I. Nuclei included for the search of the pairing interaction parameters, d and 32- 





64 Gd 


66 Dy 


68 Er 


70 Yb 


72 Hf 


74W 


N for El"-"' 
N for Jo 
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86-96 
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80-104 
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80-108 
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84-110 
90-110 


86-114 
92-114 
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Fig. 1. Root-mean-square deviations of neutron even-odd mass differences and moments of inertia, 
calculated by using the non-stretched (left), single-stretched (middle), and double-stretched 
(right) residual quadrupole-pairing interactions. Th e up per panels shows the results for Xrms 
and the lower panels for X^^, see Eqs. ( |2-12| ) and (2-lJ). They are calculated as functions of 
the two parameters d and g2. Each curve is drawn with a fixed value of d (MeV), which is 
attached near the curve, as a function of g2. 



Note that the neutron even-odd mass difference has been calculated in the same 
way as the experimental data by taking the third order difference of calculated bind- 
ing energies for even and odd A'^ nuclei, where the blocking HB calculation has been 
done for odd-mass nuclei. In the Nilsson-Strutinsky method the grid of deformation 
parameters —0.08 < £2 < 0.40 and —0.08 < £4 < 0.12 with an interval of 0.04 are 
used. The Is and // parameters of the Nilsson potential are taken from Ref. 31). We 
have assumed the axial symmetry in the calculation of this subsection since only the 
ground state properties are examined. Experimental binding energies are taken from 
the 1993 Atomic Mass Evaluation. E3) As for the experimental moment of inertia, 
the Harris parameters, and Ji, are calculated from the observed excitation ener- 
gies of the 2"*" and 4"*" states belonging to the grouncLstate band, experimental data 
being taken from Ref. 33) and the ENSDF database. If the value of i7o calculated 
in this way becomes negative or J'l greater than 1000 /i^/MeV^ (this happens for 
near spherical nuclei), then j7o is evaluated by iinly using the 2"*" energy, i.e. by 
3/E2+. The Thouless-Valtion moment of inertia, which includes the effect of the 
K = 1 component of the residual quadrupole-pairing interaction, is employed as the 
calculated moment of inertia. Here, again, the matrix elements between states with 
^-^osc = i2 are neglected for simplicity in the same way as in the step of diagonal- 
ization of the mean-field Hamiltonian. The contributions of them are rather small 
for the calculation of moment of inertia, since the Z\A^osc = i2 matrix elements of 
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the angular momentum operator are smaller than the ANqsc = ones by a factor 
~ €2 and the energy denominators are larger. We have checked that those effects are 
less than 5 % for the Thouless-Valatin moment of inertia in well deformed nuclei. 

In Fig. 1 we show root-mean-square deviations of the result of calculation for 
neutron even-odd mass differences and moments of inertia. We have found that the 
behaviors of these two quantities, Xrms{Ei^ and Xrms(i7o), as functions of 52 with 
fixed d are opposite, and so the mean value 

become almost constant, especially for the case of the non-stretched quadrupole- 
pairing. Therefore, we also display the results for the maximum among the two. 



XZs = max{Xrn,s(4"°^), ^rms(Jo)}. (2-13) 



As is clear from Fig. 1, the best fit is obtained for the double-stretched quadrupole- 
pairing interaction with d = 14 (MeV) and (72 = 30. It should be mentioned that the 
value of g2 is close to the one 52 = 287r/3 in Ref. [2^ ), where it is derived from the 
multipole decomposition of the (5-interaction and this argument is equally applicable 
if the double-stretched coordinate is used in the interaction. It is interesting to notice 
that if the non-stretched or the single-stretched quadrupole-pairing interaction is 
used, then one cannot make either X^^^s or X^^ smaller than 0.2. Xrms in the non- 
stretched case is rather flat as a function of 52 and the minimum occurs at d = 12 
(MeV) and 52 = (no quadrupole-pairing). X^^ in the non-stretched case takes 
the minimum at small quadrupole-pairing, d = 12 (MeV) and 52 = 2. Both Xj-ms 
and are flat as a function of §2 also in the single-stretched case, and take the 

minimum at d = 12 (MeV) and (72 = 16. In contrast, the double-stretched interaction 
gives well developed minima for both Xrms and X^^. These results clearly show 
that one has to use the double-stretched quadrupole-pairing interaction. 

One may wonder why the non- and single-stretched interactions do not essen- 
tially improve the root-mean-square deviations. The quadrupole-pairing interaction 
affects eI^'"^ and Jo in two ways: One is the static (mean-field) effect through the 
change of static pairing potential ( |2-5D , and the other is a dynamical effect (higher 
order than the mean-field approximation) and typically appears as the Migdal term 
in the Thouless-Valatin moment of inertia (c.f. Eq. ( ^-73 ) and ( |3-74 )). The former 
effect can be estimated by the averaged pairing gap, 

i i K i i 

where the summation is taken over the Nilsson basis state i's included in the pairing 
model space. Stronger quadrupole-pairing interaction results in larger A, which leads 
to the increase of even-odd mass difference on one hand and the reduction of moment 
of inertia on the other hand. The Migdal term coming from the K = 1 component 
of the quadrupole-pairing interaction makes the moment of inertia larger when the 
force strength is increased. Therefore, the moment of inertia either increases or 
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Fig. 2. Pairing gaps (upper panels) and moments of inertia (lower panels), calculated by using the 
non-stretched (left), single-stretched (middle), and double-stretched (right) residual quadrupole- 
pairi ng in teractions. Average pairing gaps A and monopole-pairing gaps Zioo (MeV), see 
Eq. (2-14), are displayed by solid and dashed curves, respectively, in t he u pper panels, while 
Thouless-Valatin and Belyaev moments of inertia {Ti^ /MeV) (c.f. Eq. (3-73)) are displayed as 
solid and dashed curves, respectively, in the lower panels. They are calculated as functions of the 
two parameters d and g^- Each curve is drawn with a fixed value of d (MeV), which is attached 
near the curve and changed by step of 1 MeV, as a function of g2. The calculation has been done 
for a typical deformed nucleus, '^^^Yb, with deformation parameters (e2, £4) ~ (0.2570, 0.0162). 



decreases as a function of force strength, depending on which effect is stronger. 
In Fig. 2, we have shown the energy gap and the moment of inertia for a typical 
rare-earth deformed nuclei ^^^Yb as functions of the two parameters d and 52 in 
parallel with Fig. 1. One can see that the average as well as monopole-pairing gaps 
increase rapidly as functions of the quadrupole-pairing strength if the non-stretched 
interaction is used. This static effect is so strong that the Thouless-Valatin moment 
of inertia decreases. In the case of the single-stretched case, similar trend is observed 
for the pairing gap, though it is not so dramatic as in the case of non-stretched 
interaction. The static effect almost cancels out the dynamical effect and then the 
Thouless-Valatin moment of inertia stays almost constant against (72 in this case. 
On the other hand, if one uses the double-stretched interaction, the pairing gap 
stays almost constant as a function of §2- This is because {P2K) ~ holds in a 
very good approximation, which is in parallel with the fact that the quadrupole 
equilibrium shape satisfies the selfconsistent condition, {Q2K) = 0, for the double- 
stretched quadupole operator. Thus the effect of the double-stretched quadrupole- 
pairing interaction plays a similar role as the particle-hole interaction channel; it 
acts as a residual interaction and does not contribute to the static mean-field. 
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2.4. Results of calculation 

It has been found in the previous subsection that the double-stretched form 
of the quadrupole-pairing interaction with parameters d = 14 MeV and g2 = 30 
gives best fitting for the even-odd mass differences and the moments of inertia in 

the rare-earth region. Resulting root-mean-square deviations are Xr^is{E^''^\ J'q) = 
(0.115, 0.136). If one uses {d,g2) = (13,28) or (12,20), as examples, those quantities 
become X,^s{Ei'''°\ Jo) = (0.154, 0.127) or (0.235, 0.121), respectively. Therefore, 
making the two quantities smaller is complementary as discussed in §2.3. 




90 100 

N 



110 120 



Fig. 3. Comparison of calculated even-odd mass differences (left panels, in MeV) and moments of 
inertia (right panels, in h'^ /MeV) with experimental data for rmclei in the rare-earth region. 
Experimental data are displayed in the upper panels while the calculated ones in the lower 
panels. Isotopes with Z = 62 — 78 are connected by solid {Z = mod 4) or dashed {Z = 2 
mod 4) curves as functions of neutron number N. The double-stretched quadrupole-pairing 
interaction is used with parameters d = 14 MeV and g2 = 30. 



We compare the results of calculation with experimental data in Fig. 3 as func- 
tions of neutron number. In this calculation the results of Sm {Z = 62), Os {Z = 76) 
and Pt {Z = 78) isotopes are also included, which are not taken into account in the 
fitting procedure. As is clear from the figure, both even-odd mass differences and 
moments of inertia are not well reproduced in heavy Os and Pt isotopes; especially 
even-odd mass differences are underestimated by about 20%, and moments of inertia 
overestimated by about up to 50% in Pt nuclei with N > 100. In these nuclei, low- 
lying spectra suggest that they are 7-unstable, and therefore correlations in the 7 
degrees of freedom are expected to play an important role. Except for these nuclei, 
the overall agreements have been achieved, particularly for deformed nuclei with 
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Fig. 4. Even-odd mass differences (left panels) and moments of inertia (right panels) for Z = 64— 74 
isotopes, calculated by using the single- and non-stretched quadrupole-pairing interactions. The 
panels from top to bottom show the results of the single-stretched cases with parameters (d = 11 
MeV, g2 = 18) and (d = 12 MeV, g2 = 16), and of the non-stretched cases with parameters 
(d = 11 MeV, g2 = 5) and (d = 12 MeV, 52 = 0), respectively. 
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Fig. 5. Same as Fig. 3, but for nuclei in the actinide region. Isotopes with = 88 — 100 are 
connected by solid (Z = mod 4) or dashed {Z = 2 mod 4) curves as functions of neutron 
number N . 



N K 90— 100. It is, however, noted that some features seen in experimental data are 
not reproduced in the calculation: (1) The maximum at = 90 and the minimum at 
A*" = 106 or 110 in the even-odd mass difference are shifted to = 92 and N = 98, 
respectively. This is because details of the neutron single-particle level spacings in 
the present Nilsson potential are slightly inadequate. (2) The proton number de- 
pendences of both the even-odd mass difference and the moment of inertia are too 
weak: curves of both quantities bunch more strongly in the calculation. This trend 
is clearer in light Z nuclei, Z < 68, for example, Gd or Dy; the even-odd mass dif- 
ference in these isotopes decreases more slowly as a function of neutron number in 
the calculation, which results in the slower increase of the moment of inertia. This 
problem suggests that some neutron-proton correlations might be necessary. 

For comparison sake, results obtained by using the quadrupole-pairing interac- 
tions of the single-stretched and the non-stretched types are displayed in Fig. 4. In 
the calculation of the single-stretched case, the values of the two parameters, d = 11 
MeV and 52 = 18, are employed, resuhing Xrms{Ei'''°\ Jq) = (0.240, 0.170), in 

one case, and the values d = 12 MeV and §2 = 16, resulting Xrms{Ei^~°\ Jo) = 
(0.192, 0.214), in another case. Comparing with the experimental data in Fig. 3, 
the decrease of even-odd mass difference with neutron number is too strong, while 
the increase of moment of inertia near A*" ~ 90 is too slow. In the calculation of 
non-stretched case, the values of the two parameters, d = 11 MeV and 52 = 5 are 
employed, resulting X,^{e'^'°\ Jq) = (0.257, 0.237), in one case, and the values 
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d = 12 MeV and 52 = 0, resulting Xrms{Er°' , Jo) = (0.265, 0.204), in another 
case. The average values of the even-odd mass difference are considerably smaller 
and those of the moment of inertia are 20—30% smaller compared to the experi- 
mental data. Note that the last case {d = 12 MeV and 52 = 0) is nothing but the 
calculation without the quadrupole-pairing interaction. The trend of weak proton 
number dependence does not change for all three forms of the quadrupole-pairing 
interaction. 

The merit of the Nilsson-Strutinsky method is that a global calculation is pos- 
sible once the mean-field potential is given. We have then performed the calculation 
for nuclei in the actinide region with the same pairing interaction and the parameters 
as in the rare-earth region, i.e. the double-stretched quadrupole-pairing with d = 14 
MeV and g'2 = 30. The result is shown in Fig. 5. Nuclei in the light actinide region 
are spherical or weakly deformed with possible octupole deformations. The exper- 
imental moments of inertia suggest that nucleus in this region begins to deform at 
N ~ 134, and gradually increases the deformation until a rather stable deformation 
is established at > 140. In the nuclei with = 88 and 90, the neutron number at 
which the deformation starts to grow is too large in the calculation, and the even- 
odd mass differences take considerably different behaviour from the experimental 
data. This disagreement possibly suggests the importance of octupole correlations. 
Except for these difficiencies, both even-odd mass differences and moments of inertia 
in heavy well-deformed nuclei are very well reproduced in the calculation. It should 
be emphasized that the parameters fixed in §2.3 for the rare-earth region are equally 
well applicable for the actinide region. 

§3. sec method for constructing diabatic rotational bands 

The sec method 0^ is a theory aiming at microscopic description of large am- 
plitude collective motions in nuclei. The rotational motion is one of the most typical 
large amplitude motions. Therefore it is natural to apply the SCC method to the 
nuclear collective rotation. In Ref. ^), this line has been put into practice for the 
first time in order to obtain the diabatic rotational bands, where the interband in- 
teraction associated with the quasiparticle alignments is eliminated. It has also been 
shown that the equation of path in the SCC method leads to the selfconsistent crank- 
ing model in the case of rotational motion. Corresponding to the uniform rotation 
about one of the principal axes of nuclear deformation, the one-dimensional rotation 
has been considered as in the usual cranking model. We keep this basic feature in 
the present work. 

More complete formulation and its application to the ground state rotational 
bands (i^-bands) in realistic nuclei have been done in Ref. |37| ), followed by further 
applications to the Stockholm bands (s-banda^Ea^ and improved calculations with 
including the quadurpole-pairing interaction. e3) In these works the basic equations 
of the SCC method have been solved in terms of the angular momentum expansion 
(/-expansion). Thus, the A and B parameters in the rotational energy expansion, 
E{I) = AI{I + 1) + B [/(/ -I- 1)]^, have been studied in detail. It is, however, well 
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known that applicability of the /-expansion is limited to relatively low-spin regions. 
This limitation is especially severe in ttie case of the s-bands: One has to take the 
starting angular momentum Iq (~ 10?iE3)) and the expansion in terms of (/ — Iq) is 
not very stable. Because of this problem comparisons with experimental data have 
not been possible for the s-bands.E3) In the present study, th^^rotational frequency 
expansion is utilized instead, according to the original work.EB) Then the diabatic 
cranking model is naturally derived. Thus, after obtaining the diabatic quasiparticle 
states, we construct the s-band as the two quasiparticle aligned band on the vacuum 
g-band at given rotation frequencies. This is precisely the method of the cranked 
shell model, which has been established as a powerful method to understand the 
high-spin rotational bands accompanying quasiparticle excitations. 



Another important difference of the present work from Refs. 38),^9|) is that 
the expansion method based on the normal modes of the random phase approxima- 
tion (RPA) is used for solving the basic equations in these references. The method is 
very convenient to investigate detailed contents of the rotation-vibration couplings, 
e.g. how each normal mode contributes to the rotational A and/or B parameters, 
as has been discussed in Refs. ^^,^8|). On the other hand, we are aiming at a sys- 
tematic study of rotational spectra of both g- and s-bands in the rare-earth region. 
Then the use of the RPA response-function matrix is more efficient for such a pur- 
pose, because it is not necessary to solve the RPA equation for all the normal modes 
explicitly. 

It has to be mentioned that the problem of nucleon number conservation, i.e. 
the pairing rotation, can be treated similarly, ej) Actually, if the SCO method is 
applied to the spatial rotational motion, the mean value of the nucleon number 
changes as the angular momentum or the rotational frequency increases. A proper 
treatment of the pairing rotations is mnuired, i.e. the coupling of the^patial and 
pairing rotations should be included. clP However, it has been foundElP that the 
efi^ect of the coupling is negligibly small for the case of the rotational motion in well 
deformed nuclei. Therefore, we simply neglect the proper treatment of the nucleon 
number in the following. 

Although it is not the purpose of this paper to review applications of the SCO 
methodto other nuclear structure phenomena, we would here like to cite a brief 
reviewE^P and some papers, in which low-frequency quadrupole vibrations are an- 
alyzed on the basis of the SCO method: anharmonic gamma vibrations, ©'El^'El* 
shape phase transitions in Srn_isotopes,@'0'El' anharmonicities of the two phonon 
states in Ru and Se isotoii^,Ej) single-particle levels and configurations in the shape 
phase transition regions, t2P and a derivation of the Bohr-Mottelson type collective 
Hamiltonian and its application to transitional Sm isotopes. ElP 

3.1. Basic formulation 

The starting point of the SCO method is the following time-dependent Hartree- 
Bogoliubov (TDHB) mean-field state 



\cl>{9,I,))=W{e,I,Mo), 



(3-1) 
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which is parametrized by the time-dependent collective variables 0{t) and Ix{t) 
through the unitary transformation W{9,Ix) from the ground (non-rotating) state 
|i;^o)- In the case of rotational motion, corresponds to the angular momentum 
about the rotating axis x, which is a conserved quantity, and 6 is the conjugate 
angle variable around the x-axis. In order to guarantee the rotational invariance, 
W{9,Ix) has to be of the form 

W{e,Ix) = e-'^'^''e'^'^^^\ (3-2) 

where Jx is the angular momentum operator about the x-axis, and G{Ix) is a one- 
body Hermite operator by which the intrinsic state is specified: 

mjx)) = e-*'-^l0i„tr(4)), |</'intr(4)) = e^^^^^^ | </)o) • (3-3) 

The generators of the unitary transformation W{9,Ix) are defined by {dW / dq)W~^ 
for q = 9 OT Ix, and they have, from Eq. ( |3-2| ), the form 

ulx Olx 

dW 1 

= Jx, (3-5) 

Oiie of the basic equations of the SCC method is the canonical variable condi- 
tions, tj) which declare that the introduced collective variables are canonical coor- 
dinate and momentum. In the present case they are given as 

{cl){e,Ix)\i0{IxmOJx)) = 0, (3-6) 

{(P{e,ix)\Jxm,ix)) = ix, (3-7) 

and form which the week canonical variable condition is derived: 

{He,Ix)\ [Jx,i0{lx)] \<P{0,Ix)) = 1. (3-8) 

The other basic equations, the canonical equations of motion for the collective vari- 
ables and the equation of path, are derived by the TDHB variational principlej*]l 

6{ci>{ejx)\{H-i-^)\cp{e,ix)) = o, (3-9) 

or by using the generators, Eqs. (3-4) and (3-5), 

^,lx)\[0, H - 9Jx + ix0{Ix)]\(l>{9,Ix)) = 0, (3-10) 



where O is an arbitrary one-body operator. Taking the generators as O and using the 
canonical variable conditions, Eqs. (3-6) — (l3-8|) , one obtains the canonical equations 
of motion: 

9 = ^=UJrot{Ix), (3-11) 

^^ = -^ = 0, (3-12) 



*•* In this subsection the unit of ?i = 1 is used. 
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with 

n{h) = {(t>{e,h)\H\(t){9,h)) = (<^intr(4)|^|0intr(4)>, (3-13) 

where the rotational invariance of the Hamiltonian, [H, Jx\ = 0, is used. Equa- 
tion (3-12) is nothing else than the angular momentum conservation, and Eq. (3-11) 
tells us that the rotational frequency is constant, i.e. the uniform rotation. Making 
use of these equations of motion, the variational principle reduces to the equation of 
path 

5(0intr(^)| H - LOrot{Ix)Jx |0intr(^)) = 0, (3-14) 

namely it leads precisely to the cranking model. The remaining task is to solve 
this equation to obtain the operator iG{Ix) under the canonical variable conditions, 
which are now rewritten as 

(0intr(4-)|C(4)|(/>intr(IxO) = 0, C(4) = ^^e'^^^^^), (3-15) 

Olx 

{(t>irAr{Ix)\Jx\(l>irAr[Ix)) = h- (3-16) 



In Ref. ^), Eqs. ( |3T4 ) — (3-16) are solved by means of the power series expansion 



method with respect to Ix, which gives the functional form of the rotational frequency 
^mt{Ix)- It is, however, well known that the convergence radius of the power series 
expansion with respect to LOmt is much larger, so that the applicability of the method 
can be enlarged. Ej) Thus, the independent variable is changed to be uj^t instead 
of Ix in the equations above. In the following, we write the rotational frequency as 
uj in place of u^ot for making the notation simpler. Now the basic equations can be 
rewritten as 

5(0i„tr('^)|-f^ - ^JxWmtri^)) = 0, (3-17) 

(0i„tr(a;)|C(a;)|0i„tr(4)) = 0, C{u:) ^ e-*^(-), (3-18) 

(0intr(t^)|-^x|'Aintr(t^)) =Ix{^)- (3-19) 

Note that the last equation is not the constraint now, but it just gives the functional 
form of the angular momentum Ix in terms of lv. The first two equations, Eqs. (3T7) 
and (3T8), are enough to get iG{Lo), which makes the calculation simpler. The 
equation of motion is transformed to the canonical relation 

^ = (3-20) 
with the total Routhian in the rotating frame 

n'iuj) = {<PinU^)\H - UjJx\<l)uM). (3-21) 
In order to show this, we note the following identity. 
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for an arbitrary u-independent one-body operator O. Then, 
dl-L' 

>intr('^)| [H - UjJx, C{uj)] l^mtrl^^)) " (<?!'intr (t^) | ^a; | </'intr (^)) (3-23) 



iG^'^Hu) = u;^{Y,g^^Hij)4a] - h.c.}. (3-25) 



which lead to Eq. (|3-2C1| ) because the first term of the right hand side vanishes due 
to the variational equation (3-17). 

The one-body operator iG{uj) generates the unitary transformation from the 
non-rotating (ground) state |(/>o), see Eq. (|3-3| ), and it is composed of the a|aj and ajtti 

terms, where o| and Oj are the creation and annihilation operators of the quasiparticle 
state i with respect to the ground state |(/)o) as a vacuum state. The solution of the 
basic equations is obtained in the form of power series expansion 

oo 

iGH = (3-24) 

n=l 

with 

It is convenient to introduce a notation for the transformed operator, which is also 
expanded in power series of w, 

oo 

b{uj) = e-iGH Oe*G(c.) ^ ^ oW(a;), (3-26) 

n=0 

for which the following formula are useful; 

oo 

e-«Oe*G = ^- [•••[0,zG]---iG], (3-27) 

n=0 ■ ^""^ ' 

n times 

and 

n=0 ^ ' 

n times 

Then the basic equations for solving iG[Lo) in the n-th order in oj are 

{U [a^ai, ^(") - |<Ao> = 0, (3-29) 

(</'o|c(")|0o) = O, (3-30) 
and the canonical relation is 

a-iy/(n+l) 

= or (n + l)W'("+i)=-u;/W, (3-31) 

where the total Routhian and the expectation value of the angular momentum are 
also expanded in power series. 



n'{u:) = Y,n'^^\ 4H = ^lW. (3-32) 



n=0 n=l 
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The lowest order solution is easily determined: The n = and 1 parts of 
Eq. (3-30) are satisfied trivially, while the n = 1 part of Eq. (3-29) is written as 



[ajau [H, iG«]]|(/.o 



uj{(po\ [ajai, 



or 



[H, «G*^^)]rpa = UjJxRPA, 



(3-33) 



(3-34) 



where the subscript [ ]rpa means that only the RPA order term is retained; e.g. 



JxRPA = and ttjUi parts of Jx- This is the RPA equation,!^ with respect 
to the ground state \4>o), for the angle operator i&rpa conjugate to the symmetry 
conserving mode J^RPA; and we obtain 



(3-35) 



where J^o is the Thouless-Valatin moment of inertia. Note that the general solution 
of Eq. ( p-33| ) contains a term iujcjJxRpa with cj being an arbitrary real constant. 
We have chosen cj = as a physical boundary condition, because Jx operator 
generates the transformation from the intrinsic to the laboratory frame and should 
be eliminated from the unitary transformation generating the intrinsic state, see 
Eq. ( |3-3D . Once the lowest order solution (n = 1) is obtained, higher order solutions 
(n > 2) can be uniquely determined by rewriting Eqs. (3-29) and (3-30) in the 
following forms; 



bo\[a,ai, [H, iG(")]]|</.o) 
i.o|[iG("),i0RPA]|<Ao) = 



Pol [ajCLi 
1 



(n - 1) Jo 



I 



with 



(n-1) 



duj 



duj 



(3-36) 
(3-37) 

(3-38) 
(3-39) 



Here B^"^^ and D^"^^ only contain iG^'^'> with m < n — 1, and diG^"^^ /div = n iG^"^^ juj 
and Eq. ( t3-35|) are used. Equation (3-36) has the same structure as Eq. ( |3-33|) 
or (3-34) and is an inhomogeneous linear equation for the amplitude g^"\ij), where 
the inhomogeneous term is determined by the lower order solutions (see §3.3 for 
details) . 

As in the case of the first order equation, if iG^"^^ is expanded in terms of the 
complete set of the RPA eigenmodes which is composed of the non-zero normal 
modes and the zero mode (J^rpAi i&KPA), tlie general solution of iG^") contains the 
term proportional to JxRpa, and it is determined by Eq. (3-37). Once the boundary 
condition for iG^^^ is chosen as above, however, the term proportional to J^rpa 
should vanish. In order to show this, one has to note that matrix elements of the 
Hamiltonian and of the angular momentum can be chosen to be real with respect to 
the quasiparticle basis (a|, Oj) in a suitable phase convention, e.g. that of Ref. |2^. 



18 



Then the matrix elements of the RPA normal mode operators and the angle operator 
i0RPA are also real, and so does the matrix elements of iG^^^ . If iG^"'^ is expanded 
in terms of the RPA eigenmodes, the imaginary part of its matrix elements arises 
only from the term proportional to J^rpa because iG^""^ is anti-Hermite while J^rpa 
is Hermite. If we assume that with m < n — 1 has no J^rpa term so that its 

matrix elements are real, then the right hand side of Eq. (3-37) vanishes, because 
Z)^") is an anti-Hermite operator with real matrix elements composed of iG^"^^ with 
m < n — 1. Therefore, iG^"^ neither contains the J^rpa term. Thus, the fact that 
the operator iG has no J^rpa term is proved by induction. The situation is exactly 
the same for the case of gauge rotation; the A^rpa term (A^ is either the neutron or 
proton number operator) also does not appear in iG. The method to solve the above 
basic equations for our case of the separable interaction ( |2T| ) will be discussed in 
detail in §3.3. 

3.2. Diabatic quasiparticle states in the rotating frame 

In the previous subsection the rotational motion based on the ground state |0o) 
is considered in terms of the SCC method. The same treatment can be done for 
one-quasiparticle states. The one-quasiparticle state is written in the most general 
form as 

\4>i-<,.M) = e*''^") E M^)4 \<Po), (3-40) 

i 

where iG{uj) as well as the amplitudes fi{uj) are determined by the TDHB variational 
principle. Generally iG{uj) for the one-quasiparicle state is not the same as that of 
the ground state rotational band because of the blocking effect. However, we neglect 
this effect and use the same iG{u)) in the present work following the idea of the 
independent quasiparticle motion in the rotating frame, cd^ Then by taking the 
variation 

^ (^l-q.p.(^)l-H' - UJJx\4'l-c^.x>i^)) 
[ (</'l-q.p.(w)|01.q.p.(w)) 

with respect to the amplitudes /j, one obtains an eigenvalue equation, 

E e[.j[u:)f,^{u) = fUuj)E'^iuj), (3-42) 
j 

with 

e^jiuj) = {(t>o\ai(^H{uj) - cjX(a;))a]|0o)- (3-43) 

Namely the excitation energy E'^{uj) and the amplitudes of the rotating quasi- 

particle state yu are obtained by diagonalizing the cranked quasiparticle Hamiltonian 
defined by 

h'iio) = one-body part of [e-'^^'^^H - ujj^) e'^^'^^] 

= E6^MaJa„ (3-44) 



(3-41) 
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where, due to the equation of path, Eq. (3-17) or (3-29), h!{oj) has no a)a^ and aa 
terms. Introducing the quasiparticle operator in the rotating frame, 

a^iu^) = e^^H M^)4 e-'^'^^l (3-45) 

i 

we can see that the one-quasiparticle state ( |3-40| ) is written as 

I'Al-q.pX'^)) =«j.('^)l</'intr(w)), (346) 

and 

h'{uj) = one-body part of {H — ojJx) 

= ^i?;Hat(a;)Q^M. (3-47) 

Namely, the quasiparticle states in the rotating frame are nothing but those given 
in the selfconsistent cranking model. Thus, if H contains residual interactions, the 
effects of change of the mean-field are automatically included in the quasiparticle 
Routhian operator (3-44) in contrast to the simple cranked shell model where the 
mean-field parameters are fixed at w = 0. 

It is crucially important to notice that the cutoff of the power series expansion 
in evaluating Eq. (3-44) results in the diabatic quasiparticle states; i.e. the positive 
and negative quasiparticle solutions do not interact with each other as functions of 



the rotational frequency. This surprising fact has been found in Ref. 36) and utilized 
in subsequent various applications to the problem of high-spin spectroscopy; see e.g. 
Ref. 52). Thus, we use 



max 



[h'{u)t^--^ = E( E (3-48) 

ij n=0 



with 

e'!:f ^ (</)o|a,(F(") -a;X("-^))a]|0o)/^", (3-49) 
as a diabatic quasiparticle Routhian operator. If we take nmax = 1 and use the 



solution ( 3-35 ), the first order Routhian operator is explicitly written as 

[^'(u;)]("<i) = h- iviJx - J^rpa), (3-50) 
with h = one-body part of H. This Hamiltonian was used to construct a diabatic 



quasiparticle basis in Ref. to study the g-s band crossing problem. We will show 
in §3.4 that the inclusion of higher order terms improves the quasiparticle Routhian 
in comparison with experimental data. 

In order to study properties of one-body observables in the rotating frame, for 
example, the aligned angular momenta of quasiparticles, an arbitrary one-body op- 



erator O has to be expressed in terms of the diabatic quasiparticle basis ( 3-45| ); 
O = e^^('^)0(w)e"*^(^) 
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fj,u 

+ ^) "m"!/ + OA-ifJ-i^; uj) a^ttf,^ (3-51) 



where the matrix elements are written as 



' f max 

OB(/xi^;^) = E/rMH/^-H('^o|ai(E o(")M)a]|0o>, (3-52) 

ij n=0 



"max ^ 

0^+(/.i.;a;) = E/*,(^)/;,(a;)(0o| [a,a,,, O^^^Hj] |0o>, (3-53) 
0^_(^i.;a;) =E/,^(a;)/,,H(0o| (E 0^"^(^))' 44^\M- (3-54) 



n=0 



It is clear from this expression that there are two origins of the w-dependence of the 
matrix elements; one is the effect of collective rotation, Eq. ( p-26 ), which is treated 



in the power series expansion in uj and truncated up to n-max) and the other comes 
from the diagonalization of the quasiparticle Routhian operator, Eq. ( |3-42| ). Our 
method to calculate the rotating quasiparticle states can be viewed as a two-step 
diagonalization; the first step is the unitary transformation e^^^^^ , which eliminates 

the dangerous terms, the a^a^ and aa terms, of the Routhian operator h' up to the 
order Umax in ^ leading to the diabatic basis, while the second step diagonalizes 
its one-body part, the a^a -terms. We shall discuss this two-step transformation in 
more detail in §4.1. In this way we can cleanly separate the effects of the collec- 
tive rotational motion on the intrinsic states of the (/-band and on the independent 
quasiparticle motion in the rotating frame. As long as the one-step diagonalization 
is performed as in the case of the usual cranking model, this separation cannot be 
achieved and the problem of the unphysical interband mixing is inevitable. 

3.3. Solution of the equation of path by means of the RPA response function 

Now we present a concrete procedure to solve the equation of path, Eq. (3-17), for 
our Hamiltonian which is composed of the Nilsson single-particle potential and the 



multi-component separable interaction (2-1). Let us rewrite our total Hamiltonian 
in the following form: 

H = h-\Y,XpQpQp, (3-55) 
where Qp are Hermite operators satisfying 

Qp = Ql, {MQp\M = ^, (3-56) 
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and |<^o) is the HB ground state of The mean-field Hamiltonian h includes the 
pairing potential and the number constraint term as well as the Nilsson Hamiltonians: 

h = %ils -EE ^LOr {Pjl + Plo) - E ^-^r, (3-57) 
T L=0,2 T 

where the nuclei under consideration are assumed to be axially symmetric at w = 0. 
Our Hamiltonians has a symmetry with respect to the 180°-rotation around the 
rotation-axis (x-axis), the quantum number of which is called signature, r = e~*"; 
therefore the operators Qp are classified according to the signature quantum num- 
bers,™ r = ±1 or a = 0, 1. Moreover, we can choose the phase conventioncSP in such 
a way that the matrix elements of the Hamiltonians H and of the angular momen- 
tum Jx are real. Then the operators Qp are further classified into two categories, i.e. 
real and imaginary operators, whose matrix elements are real and pure imaginary, 
respectively. Since expectation values of the signature r = —1 (a = 1) operators and 
of the imaginary operators vanish in the cranking model, operators with signature 
r = +1 and real matrix elements only contribute to the equation of path for the 
collective rotation. This observation is important. As shown in the end of §3.1, the 
boundary condition ( 3-35| ) for the collective rotation leads that the transformation 



operator iG{Lo) does not contain the J^rpa part in all orders. Absence of the imagi- 
nary operators guarantees that the matrix elements of iG{uj) are real and Eq. (3T8) 
is automatically satisfied: We need not use this equation anymore. 

Thus, the operators that are to be included in Eq. (3-55) in order to solve the 
basic equations for iG{oj) are 

{Qp} = -PJo+J ^20+ ) -^21+ ' -^22+ ' ^20^*322'*' (3-58) 

and correspondingly the strengths are 

{Xp] = Gl/2, Gl/2, Gl/2, Gl/2, K20, ^22, (3-59) 

where r = i/, vr distinguishes the neutron and proton operators. Here the following 
definitions are used; for the pairing operators, 

-P00+ = -Pqo + -f'oo , -Poo- = i (-Poo — -Poo 



p(±) _p(±)t , p(±) p(±) _-(p{±)\_p(±)\ .o.gnx 

^2K+ — ^2K ^ ^2K ' ^2K- — ''\^2K ^2K ) ■• W 

and for signature coupled operators, 

Qf^ = -j=L=(-.Q^K:±-Q2-K). (K>0) (3-61) 



*-* We employ the HB approximation, i.e. do not include the exchange terms of the separable 
interactions throughout this paper. 



22 



in which the superscript (it) denotes the signature r = ±1, and :0: = O — {4>q\0\4>q). 
The quasiparticle creation and annihilation operators should also be classified ac- 
cording to the signature quantum number; a\ for r = +i {a = —1/2) and a\ for 
r = —i {a = +1/2). Then the mean-field Hamiltonian h is expressed in terms of 
them as 

h = Y, (Eia\ai + Ejala^ , (3-62) 



i>0 



where X]i>o ™eans that only half of the single-particle levels has to be summed 
corresponding to the signature classification, and the quasiparticle energy at w = 
satisfies Ei = Ej. In the same way, Qp are written as 

ij>0 ij>0 

where the matrix elements satisfy, at a; = 0, qf{ji) = ■iiq^{ij) and {ij) = it q^{ij) 



for Qp with the time-reversal property being it, if the phase convention of Ref. |2S 
is used. 

Now let us consider the method to solve the equations for iG{uj). As is already 
discussed in §3.1, the solution is sought in the form of power series expansion in w, 
where the n-th order term zG^"^ is written as 



= ^" E 5^"^(u) {aH - a-ai). (3-64) 



J 

ij>0 



The n-th order equation Eq. (3-36) has the structure of an inhomogeneous linear 
equation for the amplitudes g^'^\ij), 



g('n-) \ ( 

-9^ 



where K is the RPA energy matrix 

K{ir,ki) = 



A{ij;kl) B{ij;kl) 
B*{ij;kl) A*{ij;kl) 



^ ( {M [aiai, [H, ala\] ] |</>o) {(f>o\ [ajUi, [H, ajak] ] |</'o) \ , . 
[ {M [aja], [H, 44] ] |0o) {M K4 [H, aja^] ] |<^o) J' ^ ' ^ 

and the amplitudes b^^\ij) in the inhomogeneous term are defined by 

a)a) and aa parts of ^("^ = ^ h^''\ij) {a\aX + a-jUi). (3-67) 

ij>0 

For the first order n = 1, B^^^ = ujJx^pa and Eq. ( p-65| ) determines the RPA an- 
gle operator zQrpa, as discussed in §3.1. Since the part of interaction composed of 
the imaginary operators, e.g. Pqo-, P20- and (52]"^ ^t^-' which are related to the 
symmetry recovering mode J^rpa (and N^pj^) are not included, the RPA matrix K 
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(with signature r = +1) has no zero-modes and can be inverted without any prob- 
lem. However, the dimension of the RPA matrix is not small in realistic situations, 
and therefore we invoke the merit of separable interactions; by using the response- 
function matrix for the Qp operators, the inversion of the RPA matrix is reduced to 
the inversion of the response-function matrix itself whose dimension is much smaller. 
Inserting the Hamiltonian (|3-55D into Eq. (3-36), we obtain 

{E,+Ej)g^^\ij)-Y,qj{ij)XpQi'^ =b^''\i3), (3-68) 
p 

where 

Qf^ ^ m [Qp. ^CW] |<Ao)/t^" = 2 ^ q^{ij)g'^^\ij). (3-69) 

ij>0 

in] 

Then inhomogeneous linear equations for Qp can be easily derived as 

Y.{6p„-Rp„X.)Q^^'^ =Bf\ (3-70) 

(7 

where 

_ , ^ qf{ij)q^{ij) „(n) _ ^ ^ h^''){ij)qf{ij) 

Note that Rp^ are the response functions for operators Qp and Qa at zero excitation 
energy, and nothing but the inverse energy weighted sum rule values (polarizability). 
Equation ( |3-7[1| ) is much more easily solved than Eq. ( 3-65 ) because of the huge 
reduction of dimension, and we obtain 

S^^Hu) = -Er^{T.'ifi^j)Xp[{l - i?x)-']p.4"^ (3-72) 

+ ^ pa ^ 

where the matrix notations are used for R = (Rpa) and x = i^paXp)- Apparently 
the n = 1 solution gives the Thouless-Valtin moment of inertia, 

Jo = Jty = JbcIj + Juig, JbcIj = 2 ^ [j\ ^J^ '^^ , (3-73) 

^3>0 + ^] 

and 

^Mig = E Xp[ (1 - Rxr'UBi, with ^/ ^ 2 E '^^^'^^f^'^l (3-74) 

pa ij>0 + 

where J^{ij) denote the a^a^ and aa parts of Jx, and the summation {p,cr) in 
Eq. ( p-74D runs, at a; = 0, only over p,<T = P2t+' , namely the K = I quadrupole- 
pairing component. Once the perturbative solution of iG{uj) is obtained, the quasi- 
particle energy can be calculated by diagonalizing 

h!{u:) = E [<j + ' (3-75) 

ij>0 
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and one obtains 

h'ioo) = J2 (^m(^) + E'piu^) > (3-76) 

fj.>0 

where the first and second terms in these two equations correspond to the quasipar- 
ticle states with signature r = +i (a = —1/2) and r = —i {a = +1/2), respectively. 

At the end of this subsection a few remarks are in order: First, although it is 
assumed that the starting state |0o) is the ground state at u; = 0, the formulation 
developed above can be equally well applied also when the finite frequency state at 
a; = Wo is used as a starting state; i.e. |i;^o) is determined by 6{(po\H — uJoJx\(j)o)- In 
such a case, however, the power series expansion should be performed with respect 
to {u> — ojq). In fact, the method has been applied in Ref. ^8|) to describe the s-band 
by taking the starting state as the lowest two quasineutron state at finite frequency, 
although the angular momentum expansion in (I — Iq) is used in it. Secondly, as 
can be inferred from the form of the n-th order solution (|3-72| ), the w-expansion is 
based on the perturbation with respect to the quantity uj/{Ei + Ej) (or uj/iOxiKPA), 
if the equation is solved in terms of the RPA eigenmodes). Therefore, it is expected 
that the convergence of the w-expansion becomes poor when the average value of 
the two quasiparticle energies is reduced: It is the case for the situation of weak 
pairing, or when one takes the starting state at a finite frequency where highly 
alignable two quasiparticle states have considerably smaller excitation energies. The 
difficulty in the calculation of s-band in Ref. |3^) is possibly caused by this problem. 
Thirdly, as mentioned already, the expectation value of the nucleon number is not 
conserved along the rotational band. This is because the number operator Nr does 
not commute with iG{uj); namely, there exists a coupling between the spatial and the 
pairing rotations. In order to achieve rigorous conservation of nucleon numbers, one 
has to apply the SCC method also to the pairing rotational motion, clP and combine 
it to the present formalism. In view of such a more general formulation, the energy 
in the rotating frame ( |3-21| ) calculated in the present method is actually the double 
Routhian 7i"{uj, Xr = Xqt), where Aor is the chemical potential fixed to conserve the 
number at the ground state w = 0. The w-dependence of the expectation value of 
number operator starts from the second order, and its coefficient is very small as 
will be shown in §3.4. Therefore the effect of number non-conservation along the 
rotational band is very small; this fact has been checked in Ref. 37) by explicitly 



including the coupling to the pairing rotation. Finally, this method utilizing the 
response- function matrix can be similarly applied to the case of the (r/* , 7y)-expansion 
of the SCC method for problems of collective vibration. In such full RPA 

response matrix (containing both real and imaginary operators) is necessary, and 
one has to choose one of the RPA eigenenergies, to which the solution is continued 
in the small amplitude limit, as the excitation energy of the response function. 

3.4. Application to the g- and s- bands in rare-earth nuclei 

We apply the formulation of the SCC method for the collective rotation de- 
veloped in the previous subsections to even-even deformed nuclei in the rare-earth 
region. In this calculation, the same Nilsson potential (the Is and parameters from 
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Ref. ^Tj)) is used as in fQ, but the hexadecupole deformation is not included. As in- 
vestigated in Ref. 37), the couphngs of collective rotation to the pairing vibrations 
as well as the collective surface vibrations are important. Therefore the model space 
composed of three oscillator shells, A'osc = 4 — 6 for neutrons and A'osc = 3 — 5 for 
protons, are employed and all the ANqsc = 0, ±2 matrix elements of the quadrupole 
operators are included in the calculation. In order to describe the properties of de- 
formed nuclei, the deformation parameter is one of the most important factors. The 
Nilsson-Strutinsky calculation in §^ gives slightly smaller values compared with the 
experimental data deduced from the measured B{E2, 2^ — > 0^) values. Therefore, 



we take the experimental values for the €2 parameter from Ref. 54). There exist 



however, some cases where no experimental data are available. Then we take the 
value obtained by extrapolation from available data according to the scaling of the 
result of our Nilsson-Strutinsky calculation in for example, e2(^^^Dy) used is 
e2(^^'^Dy)=^i X e2(^^^Dy)^^P/e2(^^'^Dy)^^i. The values adopted in the calculation are 
listed in Table II. 

Table II. Summary of the calculated results and comparison with experimental data for nuclei in 
the rare-earth region, Gd {Z = 64) to W (Z — 74). The deformation parameters 62 are taken 
from Ref. superscript * denotes cases where no data is available and extrapolation based 
on our calculation in §^ is employed. The Harris parameters J'o and J7i are given in unit of 
ft^/MeV and ft^/MeV^, respectively. The energy gaps A are in unit of MeV. The third order 
even-odd mass differences based of the mass table of Ref. ^) are used as experimental pairing 
gaps. 





N 


£2 


Jo 


Ji 


/jcxp 
^0 


^cxp 


/i cal 




Acxp 




Gd 


88 


0.164 


11.8 


308 


8.7 




1.157 


1.424 


1.108 


1.475 




90 


0.251 


25.6 


341 


23.1 


333 


1.270 


1.169 


1.277 


1.133 




92 


0.274 


31.5 


165 


33.4 


179 


1.222 


1.097 


1.070 


0.960 




94 


0.282 


34.2 


118 


37.6 


111 


1.152 


1.060 


0.892 


0.878 




96 


0.287 


36.0 


98 


39.7 


101 


1.073 


1.030 


0.831 


0.871 


Dy 


88 


0.205* 


17.4 


134 


9.0 




1.187 


1.261 


1.177 


1.472 




90 


0.242 


24.3 


223 


20.1 


348 


1.233 


1.138 


1.269 


1.162 




92 


0.261 


29.4 


178 


29.9 


184 


1.196 


1.073 


1.077 


1.033 




94 


0.271 


32.7 


136 


34.3 


123 


1.128 


1.033 


0.967 


0.978 




96 


0.270 


34.3 


120 


37.0 


93 


1.050 


1.013 


0.917 


0.930 




98 


0.275 


36.8 


117 


40.7 


98 


0.970 


0.984 


0.832 


0.875 


Er 


88 


0.162* 


12.2 


110 


8.7 




1.105 


1.321 


1.213 


1.396 




90 


0.204 


18.6 


112 


13.0 


281 


1.153 


1.188 


1.277 


1.244 




92 


0.245 


26.2 


154 


23.1 


196 


1.165 


1.075 


1.138 


1.137 




94 


0.258 


30.3 


130 


29.0 


133 


1.105 


1.031 


1.078 


1.091 




96 


0.269 


33.5 


104 


32.6 


93 


1.028 


0.995 


1.035 


0.987 




98 


0.272 


35.8 


108 


37.1 


105 


0.951 


0.971 


0.966 


0.877 




100 


0.271 


36.4 


103 


37.5 


57 


0.919 


0.953 


0.776 


0.857 




102 


0.268 


35.4 


76 


38.1 


59 


0.907 


0.938 


0.708 


0.797 



The residual interaction is of the form given in Eq. ( ^-iD , where the double- 
stretched form factor is taken according to the discussion in §^. However, we cannot 
use the same best values obtained in ^ for the strengths of the pairing interactions, 
since the model space and the treatment of Z\A^osc = i2 matrix elements of the 
quadrupole operators are different. Here we use Gq = 20/A MeV and Gq = 2A/A 
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Table II. continued. 





iV 


£2 




jca.1 


■^0 








/vcxp 


/icxp 


Yb 


90 


0.172* 


14.4 


83 


9.1 


221 


1.124 


1.200 


1.402 


1.253 




92 


0.197* 


18.9 


119 


16.6 


204 


1.136 


1.128 


1.168 


1.180 




94 


0.218* 


23.8 


141 


23.5 


186 


1.106 


1.070 


1.137 


1.214 




96 


0.245* 


30.1 


119 


29.0 


131 


1.024 


1.012 


1.159 


1.111 




98 


0.258 


33.6 


111 


34.0 


127 


0.950 


0.981 


1.039 


0.983 




100 


0.262 


34.9 


108 


35.5 


83 


0.915 


0.959 


0.865 


0.908 




102 


0.267 


34.5 


75 


38.0 


70 


0.889 


0.938 


0.764 


0.840 




104 


0.259 


33.4 


70 


39.1 


64 


0.862 


0.926 


0.685 


0.848 




106 


0.250 


32.3 


93 


36.4 


55 


0.847 


0.918 


0.585 


0.815 


Hf 


92 


0.163* 


14.1 


90 


12.2 


178 


1.154 


1.105 


1.219 


1.260 




94 


0.181* 


17.8 


129 


17.7 


196 


1.148 


1.057 


1.175 


1.285 




96 


0.207* 


23.6 


134 


23.5 


191 


1.083 


1.004 


1.123 


1.182 




98 


0.218* 


27.0 


122 


29.3 


194 


1.032 


0.976 


1.022 


1.062 




100 


0.227 


29.5 


116 


31.2 


131 


0.986 


0.952 


0.953 


0.988 




102 


0.235 


30.6 


94 


32.7 


110 


0.935 


0.932 


0.901 


0.915 




104 


0.245 


31.4 


74 


33.8 


88 


0.867 


0.915 


0.811 


0.864 




106 


0.227 


29.2 


99 


32.1 


65 


0.867 


0.903 


0.693 


0.824 




108 


0.227 


26.9 


100 


32.1 


40 


0.898 


0.887 


0.745 


0.856 


W 


92 


0.148* 


12.1 


70 


9.4 


159 


1.159 


1.006 


1.331 


1.295 




94 


0.161* 


14.6 


100 


13.2 


182 


1.169 


0.968 


1.201 


1.142 




96 


0.179* 


18.5 


122 


17.8 


216 


1.139 


0.928 


1.146 


1.100 




98 


0.196* 


22.6 


118 


23.4 


255 


1.082 


0.899 


1.046 


1.053 




100 


0.206* 


25.4 


110 


26.3 


171 


1.032 


0.880 


1.091 


1.023 




102 


0.211* 


26.7 


99 


27.1 


134 


0.985 


0.865 


0.931 


1.027 




104 


0.214* 


27.3 


83 


28.0 


112 


0.929 


0.850 


0.884 


1.036 




106 


0.212 


26.8 


95 


28.7 


86 


0.890 


0.833 


0.802 


0.943 




108 


0.208 


24.5 


92 


29.8 


53 


0.903 


0.817 


0.814 


0.849 




110 


0.197 


21.5 


77 


26.8 


55 


0.927 


0.805 


0.720 


0.868 




112 


0.191 


19.6 


76 


24.3 


67 


0.919 


0.794 


0.793 


0.907 



MeV for the monopole-pairing interaction, by which monopole-pairing gaps calcu- 
lated with the use of the above model space roughly reproduce the experimental 
even-odd mass differences (see Eq. ( p-10 ), and note that an extra difference of the 
constant "c" in it between neutrons and protons comes from the difference of the 
model space). As for the double-stretched quadrupole-pairing interaction, we take 
92 ~ 92 ~ 2^ (^^^ by which overall agreements are achieved for the mo- 

ments of inertia. The results are summarized in Table II. Here calculated energy 
gaps A are the monopole-pairing gaps, but they are very similar to the average 
pairing gaps (2-14) because the double-stretched quadrupole-pairing interaction is 
used. The isoscalar (double-stretched) quadrupole interaction does not contribute 
to the Thouless-Valatin moment of inertia Jq, but affects the higher order Harris 
parameter J7i. We do not fit the strengths for each nucleus, but use K2k = 1.45 kI^^^ 
(see Eq. ( |2-3D ), which gives, on an average, about 1 MeV for the excitation energy of 
7-vibrations in the above model space. We believe that this choice is more suitable 
to understand the systematic behaviors of the result of calculation for nuclei in the 
rare-earth region. 
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One of the most important output quantities is the rotational energy parameters, 
i.e. the Harris parameters, in our formalism of the w-expansion. Up to the third 
order, 

4(cj) = 7 + 1/2 = Jo + (3-77) 

where / = 0, 2, 4, • • • (h) for the K = ground state bands. The results are sum- 
marized in Table II in comparison with experimental data, where the experimental 
Harris parameters Jo and Ji are extracted from the E2+ and of the ground 
state band as follows: 

_ l.buj^ - 3.5a;f _ 3.5a;i - 1.5u;2 , . 

Jo— o o , Jl — o o , l'J''oJ 

UJlUJ^ — UJ^UJ2 UJ1UJ2 — U)'Il02 

with 

LOi = E2+/2, iV2 = {Ei+ - E2+)/2. (3-79) 

If the resultant parameter becomes negative or Ji gets greater than 1000 ^^/MeV^, 
then only Jo = 3/i?2+ parameters are shown in Table II. It is seen from the table 
that two Harris parameters are nicely reproduced, especially their mass number de- 
pendence. In contrast to the Jo parameter, for which only the residual quadrupole 
pairing interaction affects, the Ji parameter are sensitive to all components of the 
residual interaction. In other words, Ji reflects the mode-mode couplings of the col- 
lective rotation to other elementary excitation modes. Therefore the SCO method 
with the present residual interaction is considered to be a powerful means to de- 
scribe the "non-adiabaticity" of nuclear collective rotations. Details of coupling 
mechanism has been investigated in Ref. ^) by decomposing the contributions from 
various RPA eigenmodes: It has been found that the couplings to the pairing vibra- 
tions and collective surface vibrations are especially important. Although the main 
contributions come from the collective modes, many RPA eigenmodes have to be 



included to reach the correct results, see also Ref. 44) for this point. The method 
of the response-function matrix described in §3.3 is very useful to include all RPA 
eigenmodes. 

Table III. Results of the a;-expansion for some observables in Er (Z = 68) isotopes. ~ 0' 2) 

are expectation values of the mass quadrupole operators. The zero-th order values of A are 
shown in Table II, and those of Q2l' 'lbyo (axially symmetric &t u = 0). Units of each 
quantity are shown in the second raw. 



iV 




ftVMeV^ 






b 


(Q2o')i 
bftVMeV^ 


b?iVMeV^ 


Er 88 


14.5 


-4.2 


-0.45 


-1.81 


2.84 


6.14 


4.58 


90 


9.9 


-5.4 


-0.72 


-1.80 


3.74 


4.67 


4.78 


92 


7.5 


-6.3 


-1.58 


-1.83 


4.71 


4.52 


5.72 


94 


3.6 


-4.3 


-2.12 


-1.50 


5.13 


2.45 


4.43 


96 


1.8 


-3.1 


-2.34 


-1.33 


5.50 


1.31 


3.70 


98 


1.0 


-3.0 


-2.83 


-1.32 


5.71 


1.10 


2.74 


100 


-2.6 


-3.4 


-2.83 


-1.35 


5.82 


0.86 


1.51 


102 


-3.5 


-3.6 


-2.34 


-1.38 


5.87 


0.79 


1.08 
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Table IV. Similar to Table 111 but the residual interactions are artificially switched off in the cal- 
culation. The results for the Harris parameters are also included. 



N 


Jo 
h'^/MeV 


Ji 


iN)i 


h'^/MeV 


(Q20*)l 

b^VMeV^ 




Er 88 


7.5 


11 


1.00 


-0.18 


0.32 


0.13 


90 


12.1 


18 


1.05 


-0.25 


0.38 


0.15 


92 


18.1 


28 


0.97 


-0.39 


0.42 


0.18 


94 


21.4 


34 


0.60 


-0.48 


0.33 


0.18 


96 


24.6 


34 


0.67 


-0.54 


0.27 


0.18 


98 


27.6 


51 


0.83 


-0.67 


0.29 


0.17 


100 


28.3 


55 


-1.01 


-0.67 


0.15 


0.15 


102 


27.5 


36 


-0.85 


-0.57 


0.16 


0.15 



Expectation values of other observable quantities are also expanded in power 
series of oj, and their coefficients give us important information about the response 
of nucleus against the collective rotation. In Table III we show examples for the 
nucleon number, monopole-pairing gaps, and mass quadrupole moments: 

(</.intr(c^)|iVr|'/'intr(c^)) = (iV,)o + {Nr)l UJ^ (3-80) 
G5(0i„tr(c^)|Po"ol'/'intr(w)) = (Z^r)o + (^r)l (3-81) 
(</'intr(a;)|Q^iVintr(^)) = (Q^i^)o + (Q^t-^)l^' {K = 0,2). (3-82) 

They are time-reversal even quantities so that the series contains up to the second 
order within the third order calculations. It should be noticed that these w-expanded 
quantities are associated with the properties of the diabatic ground state band, which 
becomes non-yrast after the g-s band-crossing. As remarked in the end of §3.3, 
{Nr)i 7^ means that the nucleon number is not conserved along the rotational 
band. However, its breakdown is rather small; even in the worst case of ^^^'Ei in 
Table III the deviation is about 1.3 at a; = 0.3 MeV, and it is less than 0.1 at uj = 0.1 
MeV in ^^^Er. It is well known that the pairing gap decreases as a function of to due 
to_the Coriolis anti-pairing effect. It is sometimes phenomenologically parametrized 
asH 

Ao(l-U—Y) uj<uj„ 
-, V 2VwcW (3.83) 

-Ao[ — w > Wc- 

Thus, our w-expansion method precisely gives the phenomenological parameter uJq = 



\J —A^jlAx < 0) in Eq. ( 3-83 ) from microscopic calculations. As shown in 
Table III, {^Ay)\ varies considerably along the isotopic chain. The {Q^k)! are related 
to the shape change at high-spin states, and tell us how soft the nucleus is against 
rapid rotation. Since nuclei studied in the present work are axially symmetric in 
their ground states, ((52o)i and (Q22)i serve as measures of softness in the /?- and 7- 
directions, respectively. As seen in Table III the isotopes get harder in both directions 
as the neutron number increases; especially, the = 88 and = 90 isotopes are 
known to undergo a shape change from the prolate collective to the oblate non- 
collective rotation scheme at very high-spin states ("band termination" B ) , while 
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heavier isotopes {N > 96) are known to be well deformed keeping prolate shape 
until the highest observed spins. These features have been well known from the 
calculations of the potential energy surface in the (e2, 7)-plane, and our results seem 
to agree with them qualitatively. In order to see the effect of the residual interactions, 
the result obtained by neglecting them, i.e. that of a simple higher order Coriolis 
coupling calculations, is shown in Table IV. Comparing it with Tables II and III, 
it is clear that the residual interactions play an important role in the w-dependence 
of observables. For example, jTi Harris parameter becomes quite small by a factor 
of about 1/2—1/3 when the residual interactions are switched off. The effect on the 
second order coefficients of the quadrupole moment is more dramatic and leads to 
about an order of magnitude reduction in soft nuclei. 




Fig. 6. Neutron quasiparticle Routhians plotted as functions of huirot (Me V) su itable for ^''^Er. 
They are obtained by diagonalizing the SCC quasiparticle Hamiltonian (3-4^) up to the first 
order (left) and third order (right) of the a;-expansion. As in the case of the usual adiabatic 
quasiparticle energy diagram, the negative energy solutions, — -E^ ~ E'fi and —E'^ = El, are 
also drawn. The solid, dotted, dashed, and dash-dotted curves denote Routhians with (tt, r) — 
(+,-i), and (-,— i), respectively. 



Now let us study the quasiparicle Routhians obtained by means of the SCC 
method. It is mentioned in §3.2 that the two step diagonalization with the truncation 
of the w-expansion up to nmaxj c.f. Eq. ( |3-48D , leads to diabatic quasiparticle states 
in the rotating frame, in which the negative and positive eigenstates do not interact 
with each other. We show in Figs. 6 and |^ calculated quasiparticle Routhians for 
neutrons and protons, respectively. It is confirmed that the diabatic quasiparticle 
states are obtained. As discussed in §3.2, the diagonalization of the quasiparticle 
Hamiltonian in the SCC method is completely equivalent to that of the selfconsistent 
cranking model, which is known to lead to the adiabatic levels, if the first step 
unitary transformation e**^^'^-* is treated non-perturbatively in full order. Then what 
is the mechanism that realizes the diabatic levels? We believe that the cutoff of the 
(j-expansion extracts the smoothly varying part of the quasiparticle Hamiltonian; 
namely, ignoring its higher order terms eliminates the cause of abrupt changes of the 
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Fig. 7. Same as Fig. 6 but for proton quasiparticles. 




Fig. 8. Same as Fig. 6 but obtained by the adiabatic cranking (left) and the third order SCC (right) 
with neglecting the residual interactions. 



microscopic internal structure by quasiparticle alignments. An analogous mechanism 
has been known for many years in the Strutinsky smoothing procedure: The 
(5-function in the microscopic level density is replaced by the Gaussian smearing 
function times the sum of the Hermite polynomials (complete set), and the lower 
order cutoff of the sum (usually 6th order is taken) gives the smoothed level density. 
It should be noted, however, that the plateau condition guarantees that the order of 
cutoff does not affect the physical results in the case of the Strutinsky method. We 
have not yet succeeded in obtaining such a condition in the present case of the cutoff 
of the w-expansion in the SCC method for the collective rotation. Therefore we have 
to decide the rimax value by comparison of the calculated results with experimental 
data. We mainly take nmax = 3 in the following; determination of the optimal choice 
of ^max remains as a future problem. 
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Fig. 9. Neutron quasiparticle Routhians for the A^'osc = 6 (ii3/2) orbits suitable for Er. The 
left upper, right upper, left lower, and right lower panels denote the results of the adiabatic 
cranking, the SCC up to the first order, 3rd order, and 5th order, respectively. The solid and 
dotted curves denote Routhians with r = +i {a = —1/2) and r = — z (a = +1/2), respectively. 
The effect of residual interaction is completely neglected in this calculation. 



In Figs. 6 and|^the results obtained by truncating up to the first order (nmax = 1) 
and the third order (n-max = 3) are compared. It is clear that the higher order terms 
considerably modify the quasiparticle energy diagrams. Especially, the alignments 
of the lowest pair of quasiparticles are reduced for neutrons (low K states of the 
z^3/2-orbitals), while they are increased for protons (medium K states of the /iii/2- 
orbitals). Thus, the higher order effects depend strongly on the nature of orbitals. 
It should be stressed that the effects of the residual interaction, i.e. changes of the 
mean-field against the collective rotation, are contained in the quasiparticle diagrams 
presented in these fi gur es. In this sense, they are different from the spectra of the 
cranked shell model, c2P where the mean- field is fixed at a; = 0. In Fig. 8 are dis- 
played the usual adiabatic quasineutron Routhians and the third order SCC Routhi- 
ans, in both of which the residual interactions are neglected completely. Again, 
by comparing Fig. 8 with Fig. 6, it is seen that the effect of residual interactions 
considerably changes the quasiparticle states. In relation to the choice of nmaxj we 
compare in Fig. 9 the Routhians obtained by changing the cutoff order n-max = 1,3,5. 
In this figure, the usual non-selfconsistent adiabatic Routhians are also displayed, 
and for comparison's sake, the residual interactions are completely neglected in all 
cases. Moreover, the rotational frequency is extended to unrealistically large values 
in order to see the asymptotic behaviors of the Routhian. Comparing the adiabatic 
Routhians with those of the SCC method, positive and negative energy solutions 
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cross irrespective of the strength of level-repulsion. Although the adiabatic levels 
change their characters abruptly at the crossing, if their average behaviors are com- 
pared to the calculated ones, the third order results (nmax = 3) agree best with the 
adiabatic levels. The first order results, for example, give the alignments (the slopes 
of Routhians) too large. On the other hand, the divergent behaviors are clearly seen 
at about uj > 0.8 MeV in the fifth order results. The inclusion of the effect of the 
residual interactions makes this convergence radius in uj even smaller. 




Fig. 10. Comparison of the third order SCC method calculations for the diabatic g- and s-bands 
with experimental data. The angular momenta (J^) = I + 1/2 (h) are displayed versus the 
rotational frequency JuOrot (MeV) for nuclei in the rare-earth region, Gd {Z = 64) toW {Z — 74) 
isotopes. Filled circles denote experimental data smoothly extended from the ground state. Data 
for excited bands are also included as filled squares when available, which are, in most cases, 
identified as s-bands. 

Finally, we would like to discuss the results of application of the present for- 
malism to the g- and s-bands, which are observed systematically and compose the 



yrast lines of even-even nuclei. Although we can compare the Routhians (3-21), or 
equivalently the rotational energy ( |3T3| ), it is known that the relation versus oj 
gives a more stringent test. Therefore we compare the calculated Ix — uJ relation 
with the experimental one in Fig. 10 for even-even nuclei in the rare-earth region, in 
which the band crossings are identified along the yrast sequences. In this calculation 
the Ix{^) of the (7-band is given by Eq. ( [j-77| ) with calculated values of the Harris 
parameters (see Table II). As for the Ix{^) of the s-band, we calculate it on the 
simplest assumption of the independent quasiparticle motions in the rotating frame, 
which is the same as that of the cranked shell model: 

\(t>s{i^)) = aI(a;)4(a;)|(Aintr(u;)>, \(t>a{u:)) = |<^intr(c^)>, (3-84) 
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Fig. 10. continued. 



where ai(a;) and Q:j(a;) are the lowest r = +i and r = —i quasineutron creation 
operators in the rotating frame. Then, the Ixi^jj) of the s-band is the sum of /a;(w) 
of the 5-band and the aUgned angular momenta of two quasineutrons, which are 
calculated according to Eqs. (3-51)— (3-54), 

^ =h{io)+h{uj)+ii{uj), (h{io)) ^ =Ix{oj), (3-85) 

V / s-band V / g-band 

or by using the canonical relation between the Routhian and the aligned angular 
momentum, the alignments and can be calculated as usual: 

Since our quasiparticle Routhians behave diabatically as functions of the rotational 
frequency, the resultant g- and s-bands are also non- interacting bands; we have to 
mix them at the same angular momentum to obtain the interacting bands corre- 
sponding to the observed bands. Such a band mixing calculation is straightforward 
in our formalism if the interband g-s interaction is provided. However, it is a very 
difficult task as long as the usual adiabatic cranking model is used. In the present 
stage we are not able to estimate the g-s interband interaction theoretically. There- 
fore, we do not attempt to perform such band-mixing calculations in the present 
paper (but see §4.2). 

Looking into the results displayed in Fig. 10, one see that our diabatic formalism 
of collective rotation based on the SCC method is quite successful. The overall 
agreements are surprisingly good, considering the fact that we have only used a 
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Fig. 10. continued. 



global parametrization of the strengths of the residual interaction: 

= 20/ A, = 24/ A (MeV), (3-87) 
g!^ = gl = 24, (3-88) 
/^2x = l•45/^f^^ (3-89) 

for the model space of three A/'osc-shells (4 — 6 for neutrons and 3 — 5 for protons). The 
agreements of the calculated i^-bands come from the fact that the Harris parameters 
(Table II) are nicely reproduced in the calculation. Further agreements of the s- 
bands are not trivial, and tell us that we have obtained reliable diabatic quasiparticle 
spectra (Figs. 6 and ^). It is known that, if the Ix — uj relations of s-bands are 
parametrized in the form, Ix = i + Jo uj + Ji uj'^ , the Ji Harris parameters of s- 
bands are systematically smaller than those of g'-bands. This feature is quite well 
reproduced in the calculations, as is clearly seen in Fig. 10, and the reason is that the 
value of the aligned angular momentum of two quasineutrons decreases as a function 
of (jj. The suitable decrease is obtainable only if the residual interactions are included 
and the diabatic quasiparticle Routhians are evaluated up to the third order. 

§4. Diabatic quasiparticle basis and the interband interaction between 

the g- and s-bands 

The formulation of the previous section gives a consistent perturbative solution, 
with respect to the rotational frequency, of the basic equations of the SCC method 
for collective rotation. However, it has a problem as a method to construct the dia- 
batic quasiparticle basis: The wave functions of the diabatic levels are orthonormal 
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Fig. 10. continued. 



only within the order of cutoff (nmax) of the w-expansion. In the previous section 
only the independent quasiparticle states, i.e. one-quasiparticle states or the g- and 
s-bands, are considered and this problem does not show up. The quasiparticle states 
have another important role that they are used as a basis of complete set for a more 
sophisticated many body technique beyond the mean-field approximation; for exam- 
ple, the study of collective vibrations at high spin in terms of the RPA method in 
the rotating frame. i^'il^'ij'il^'El''® In such an application it is crucial that the 
diabatic quasiparticle basis satisfies the orthonormal property. We present in this 
section a possible method to construct the diabatic basis satisfying the orthonormal- 
ity condition. 

Another remaining problem which is not touched in the previous section is how 
to theoretically evaluate the interband interaction between the ground state band 
and the two quasineutron aligned band. Since we do not have satisfactory answer 
yet to this problem, we only present a scope for possible solutions at the end of this 
section. 

4.1. Construction of diabatic quasiparticle basis in the SCC method 

Although the basic idea is general, we restrict ourselves to the case of collec- 
tive rotation and use the good signature representation with real phase convention, 
introduced in §3.3, for the matrix elements of the Hamiltonian H and of the an- 
gular momentum J^- First let us recall that the diabatic quasiparticle basis in the 
rotating frame is obtained by the two step unitary transformation ( p-45| ). The first 
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transformation by e*'^^'^-' can be represented as foUows, 



giG(c^)^tg-iG(<^) 



j>0 
j>0 



j>0 



COS yg^g 



j>0 



9 



sin Vff^ 
T sin y/gg^ 



a-, 



\/g<F 



(4-1) 
(4-2) 



with real matrix elements gij{uj) = X]n>i see Eq. (3-64), where denotes 

the transpose of g. Thus, by using an obvious matrix notation, the transformation 
to the rotating quasiparticle operator from the w = quasiparticle operator is given 
as 



a 
at 



f 

r 



cos 



-9 



sin Vff^ \ 



T sin y/gg^ 
9 cos y/g^g 



V ^99 



(4-3) 
(4-4) 



where the real matrix elements fifj,{uj) and fi^{uj) are the amplitudes that diagonalize 
the quasiparticle Hamiltonian in the rotating frame, see Eq. ( 3-42| ), for signature 
r = +i and —i, respectively. The cutoff of the cj-expansion means that the generator 
iG{u!), i.e. the matrix g, is solved up to the n = rimax order. 



(4-5) 



n=l 



and at the same time the transformation matrix G{uj) itself is treated perturbatively 



1 _ uj^gi^)gm + . . . 



-u;g 



(1)T 



+ 



1 _ c^2^(l)T^(l) + . . . 



(4-6) 



while the other one, J'{uj), is treated non-perturbatively by the diagonalization proce- 
dure. The origin of difficulty arising when the diabatic basis is utilized as a complete 
set lies in this treatment of G{io), because the orthogonality of the matrix G{uj) is 
broken in higher-orders. 

Now the solution to this problem is apparent: The generator matrix g{io) is 



solved perturbatively like in Eq.(4-5), but the transformation matrix G{uj) has to 
be treated non-perturbatively as in Eq. (4-3). In order to realize this treatment we 
introduce new orthogonal matrices, D and D, which diagonalize gg^ and g^g within 
the signature r = +i and —i states, respectively, 



^{99^)ijDjk = Dik9l, ^[g^g)ijDjk = Dik9l, 



(4-7) 
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where we have used the fact that the matrices gg^ and g^ g have common eigenvalues, 
which are non-negative, and then we have 

_ / D(cose)Z)T gD{sin9/e)D"^ \ 
^^'^'-\-g'^D{sme/e)D'^ ^(cos0)^T j- ^4-8) 

Here (cos0) and (sin 0/0) denote diagonal matrices, whose matrix elements are 
6ij cos 6i and 6ij sin 6i/9i, respectively. The physical meaning is that the orthogo- 
nal matrices D and D are transformation matrices from the quasiparticle operators 
(aj,aj) and (at,aj) at a; = to their canonical bases, which diagonalize the density 
matrices p and p with respect to the rotational HB state {(pmtri'^)) , respectively; 



Pij = {4>mtr{i^)\alaj\(j)inti{uj)) = COS \f gg'^ 



Pij = (0mtr('^)|a|aj|0mtr(t*^)) = COS \f^g ... (4-9) 

Thus the method to construct the rotating quasiparticle basis is summarized as 
follows. First, solve the basic equation of the SCC method and obtain the generator 
matrix g{oj) up to the Timax order as in Eq. ( |4-5D . At the same time, diagonalize 
the quasiparticle Hamiltonian and obtain the eigenstates as in Eq. (|34^ ) for both 
signatures r = ibz. Secondly, diagonalize the density matrices (4-9), or equivalently 
Eq. (|^, and obt ain the orthogonal matrices D and D of the canonical bases. 



Finally, by using these matrices D and D calculate the transformation matrix Q{uj) 
as in Eq. (4-8), and then the basis transformation is determined by Eq. (4-4). 

It is instructive to consider a concrete case of the cranked shell model; i.e. the 
effect of residual interactions or the selfconsistency of mean-field is neglected aX uo > 
0. The quasiparticle basis is obtained by diagonalizing the generalized Hamiltonian 
matrix: 

/ /iNiis - \(uv\^(uv\(e' \ 

1^ -A -[h^,,^ + u^j^) jyv U ) \V U ) \ Q -E' )^ 

(4-10) 

where /inus ^-nd jx denote matrices with respect to the Nilsson (or the harmonic 
oscillator) basis at w = 0, and (?7, V) and (C/, V) are coefficients of the generalized 
Bogoliubov transformations from the Nilsson nucleon operators (c|,Cj) and (ci,q) 
(in the good signature representation), 

«M = E(^iM4 + Vif^ci), 4 = Y.^Ui^4 + ^ifCi), (4-11) 

i>0 i>0 



or in the matrix notation 




In contrast, the transformation U is decomposed into three steps in our construction 
method of the diabatic quasiparticle basis; (i) the Bogoliubov transformation IAq 
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between the nucleon (c, c^) and the quasiparticle (a, a^) at u = 0, 



4 



7t 



u 



V 

u 



(4-13) 



0, (they are diagonal, e.g. 



UiSij, if only the monopole-pairing interaction is included), (ii) the transfor- 



where u and v are the matrices of transformation at uj 

mation matrix g{uj) in Eq. (4-4), generated by e'^^^^\ and (iii) the diagonalization 
step of the rotating quasiparticle Hamiltonian T{oj) in Eq. (4-4), see also Eq. ( ^-421) , 
namely 

U{ojf^^ =UQg{oj)T{oj). (4-14) 

Here both Q{lo) and J-{oo) depend on the order of cutoff nmax in solving the generator 
iG{ijj) by the w-expansion method, but they themselves have to be calculated non- 



perturbatively, especially for Q by Eq. (4-8) with (4-7). As noticed in the end of 
§3.3, we can apply the SCC method starting from the finite frequency ojq. In such a 
case IAq is the transformation at u; = wqj and Q and T are obtained by expansions 
in terms of {u — ujq); thus, 



U{uj) 



see 



Uq{ujq) Q{u} — ujq)J^{uj — ujq) if started at u; = ujq. 



(4-15) 



It should be stressed that the transformation ( [4-14 ) only approximately diagonalize 
the Hamiltonian in Eq. ( ^-101 ) within the nmax order in the sense of w-expansion. 
Namely, some parts of the Hamiltonian corresponding to the terms higher order than 
'^max are neglected, and this is exactly the reason why we can obtain the diabatic 
basis, whose negative and positive solutions are non-interacting. 

In the case where the effect of residual interactions is neglected, i.e. correspond- 
ing to the higher order cranking, we can easily solve the basic equations of the SCC 
method. It is useful to present the solution for practical purposes; for example 
for the construction of the diabatic quasiparticle basis for the cranked shell model 
calculations. The solutions for up to the third order are given as follows: 



9^'Kij) 



1 



Ei + E- 



(4-16) 
(4-17) 



Ei + E-j 



(4V)+5 



(2) J; 



1 



+^(4V)^5« +25«j,^v^) 



and the solutions for the rotating quasiparticle Hamiltonian ( |3-48| ) — ( 3-4g| ): 



'(0) 
'(2) 



• E- 
tB 



^ijEi, 



l(4V)T + ,«j,^T; 



^(2) 



(4-18) 

(4-19) 
(4-20) 
(4-21) 
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elf = liJ^g^'^^ + 9^'^4^h, # = \iJf9^'' + 9^'^^4h, (4-22) 

where the quasiparticle energies at the starting frequency are given in Eq. ( |3-62D , 
and the matrix elements of Jx at the starting frequency are given as in Eq. (3-63) 
with Qp replaced by Jx- If the starting frequency is a; = 0, then E'j = Ei, and 
the matrix elements of Jx satisfy the relations, Jx'^ = —J^-, Jx — ~Jxi 
J^^ = J^ . The transformation G{io) is calculated from Eqs. (4-16) — (4-18), and 
J^{ui) from Eqs. (4-19) — (4-22). It should be mentioned that the selfconsistent mean- 
field calculation is in principle possible in combination with the diabatic basis pre- 
scription presented above. 

4.2. Estimate of the g-s interaction 

Once the diabatic g- and s-bands states (|3-84 ) are obtained as functions of 



one can immediately construct them as functions of angular momentum /, because 
the Ix —ijJ relation has no singularity, as shown in Fig. 10, and can easily be inverted: 

\ct>g{l)) = = \U^S{I))), (4-23) 

where u!g{I) and ujs{l) are the inverted relations of ( ^-SSf ) with Ix = I + 1/2. Physi- 
cally, one has to consider the coupling problem between them at a fixed spin value 
/. It is, however, a difficult problem because one has to calculate, for example, a 
matrix element like {(j)siI)\H\(j)g{I)) , which is an overlap between two different HB 
states; they are not orthogonal to each other due to the difference of the frequencies 
UJg{I) and uJs{I)- Although such a calculation is possible by using the Onishi for- 
mula for the overlap of general HB states, E^) it would damage the simple picture 
of quasiparticle motions in the rotating frame, and is out of scope of the present 
investigation. 

Here we assume that the wave functions varies smoothly along the diabatic rota- 
tional bands as functions of spin / or frequency co, so that the interband interaction 
between the g- and s-bands can be evaluated at the common frequency by 

Vg-sil) = {U^gs{l))\H\<Pg{Ugs{l))), (4-24) 

where ujgs is defined by an average of ujg and ujg, 

, + (4.25) 

We note that this quantity corresponds, in a good approximation, to the crossing 
frequency u)^'^ at the crossing angular momentum 

a;,,(Ir)^<^ (4-26) 

where a;^"'* is defined as a frequency at which the lowest diabatic two quasiparticle 
energy vanishes, E[{uj) + E'^{u) = 0. Using the fact that \4>s{(jj)) is the two quasi- 
particle excited state on \(f>g{u!)) (see Eq. ( ^-841) ), the interaction can be rewritten 
as 

Vg.s{I) = u;gs{I){M^gs{mJx\(pgMl))), (4-27) 
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because of the variational principle (3-17). Applying the idea of w-expansion and 
taking up to the lowest order, we have, at the crossing angular momentum I^''^, 

v,-s{ir) ^ E Mi-^nM-^ru^iij), (4-28) 

ij>0 

where fii{io) and fjiiuj) are the amplitudes of the diabatic quasiparticle diagonal- 
ization ( |3-42| ) for the lowest r = ±i quasineutrons, and should be calculated non- 
perturbatively with respect to uj. 



In Fig. 11 (right panel), we show the result evaluated by using Eq. (4-28) for 
a simple single-j shell model (^13/2) with a constant monopole-pairing gap and no 
residual interactions, in which the single-particle energies are given by 

e, = ^^!I±zAl±}l (m, = l/2,...,j), (4-29) 
JU + 1) 

with a parameter k describing the nuclear deformation. In this figure other quan- 
tities, the alignment of the lowest two quasiparticle state, the number expectation 
value, and the crossing frequency are also shown as functions of the chemical poten- 
tial. These quantities can also be evaluated in terms of the usual adiabatic cranking 
model, and they are also displayed in the left panel. Note that in the adiabatic 
cranking model the crossing frequency is defined as a frequency at which the adi- 
abatic two quasiparticle energy E'^'^'^\io) + E'I^^'^\uj) becomes the minimum, and 
the interband interaction is identified as the half of its minimum value. 0* As is 
well known, EP the g-s interaction oscillates as a function of the chemical potential, 
and both the absolute values and the oscillating behavior of the result of calcula- 
tion roughly agree with the experimental findings. Comparing two calculations, the 
interband interaction ( [4-28 ) seems to give a possible microscopic estimate based on 



the diabatic description of the g- and s-bands. We would like to stress, however, 
that its derivation is not very sound. It is an important future problem to derive the 
coupling matrix element on a more sound ground. 



§5. Concluding remarks 



In this paper, we have formulated the SCC method for the nuclear collective 
rotation. By using the rotational frequency expansion rather than the angular mo- 
mentum expansion, we have applied it to the description of the g- and s-bands 
successfully. The systematic calculation gives surprisingly good agreements with 
experimental data for both rotational bands. It has been demonstrated that the 
resultant quasiparticle states develop diabatically as functions of the rotational fre- 
quency; i.e. the negative and positive energy levels do not interact with each other. 
Although the formulation is mathematically equivalent to the selfconsistent cranking 
model, the cutoff of the w-expansion results in the diabatic levels and its mechanism 
is also discussed. The perturbative w-expansion is, however, inadequate to use the 
resultant quasiparticle basis states as a complete set. We have then presented a 
method to construct the diabatic quasiparticle basis set, which rigorously satisfies 
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Fig. 11. The g-a interband interaction (solid), the crossing frequency h,u)gs (dash-dotted) in MeV, 
the alignment i (dotted), and the expectation value of number operator n (dashed), plotted as 
functions of the chemical potential A in MeV, for the 113/2 single- j shell model without residual 
interactions. The result of the usual adiabatic cranked shell model is displayed in the left panel, 
while that of the diabatic SCC 1st order calculation in the right panel. Here the alignment i 
and the number n is scaled by their maximum values, imax = 12?i and nmax = 14. The energy 
unit is chosen such that the splitting of the ii3/2-shell roughly reproduces that of a typical well 
deformed rare-earth nucleus; i.e. n — 2.5 MeV in Eq. (1-29), and the constant A = 1.0 MeV is 
used. 



the orthonormality condition and can be safely used for the next step calculation, 
e.g. the RPA formalism for collective vibrations at high-spin. 

In order to obtain a good overall description of the rotational band for nuclei 
in the rare-earth region, we have investigated the best possible form of residual 
quadrupole-pairing interactions. It is found that the double-stretched form factor 
is essential for reproducing the even-odd mass difference and the moment of inertia 
simultaneously. 

Since the calculated g- and s-bands in our formulation are diabatic rotational 
bands, the interband interaction between them should be taken into account for 
their complete descriptions. As in any other mean-field model, however, the wave 
function obtained in our formalism is a wave packet with respect to the angular 
momentum variable. Therefore, it is not apparent how to evaluate the interband 
interaction from microscopic point of view. We have presented a possible estimate of 
the interaction, which leads to a value similar to that estimated by the level repulsion 
in the adiabatic cranking model. Further investigations are still necessary to give a 
definite conclusion to this problem. 
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